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Abstract 

We propose a scheme allowing to observe the evolution of a quantum system in the semiclassical 
regime along the paths generated by the propagator. The scheme relies on performing consecutive 
weak measurements of the position. We show how "weak trajectories" can be extracted from the 
pointers of a series of measurement devices having weakly interacted with the system. The prop- 
p I erties of these "weak trajectories" are investigated and illustrated in the case of a time-dependent 

model system. 
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In classical physics, the evolution of a physical system is given in terms of trajectories. 
Instead quantum mechanics forbids a fundamental description based on trajectories. Nev- 
ertheless the Feynman path integral approach gives a sum over paths formulation of the 
evolution of a quantum system, and when the actions are large relative to ^ - the semi- 
classical regime-, the wavefunction evolves essentially along classical paths, those of the 
corresponding classical system ij]. Of course, this does not mean that a quantum object 
is a localized particle moving on a definite path. But trajectories may remain significant 
in quantum systems: the large scale properties, experimentally observed in many systems 
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ll, |2| , display the signatures of the underlying classical dynamics. 

In this work, we aim to go further by proposing a scheme allowing to observe the evolution 
of a quantum system in the semiclassical regime along the trajectories of the corresponding 
classical system. The scheme relies on performing consecutive weak measurements (WM). 
WM jsl, y| are characterized by a very weak coupling between the system and the mea- 
surement apparatus. Thus measuring weakly an observable A results in leaving the former 
essentially unperturbed while the latter picks up on average a limited amount of information 
encapsulated in the weak value 

/i\ -MAM). m 

\ip) is the initial ('preselected') state and (x| is the final ('postselected') state obtained 
by performing a standard strong measurement after having measured A weakly. WM are 
receiving increased attention, either as a technique for signal amplification jsl or as a tool to 
investigate fundamental problems, from a theoretical standpoint but also experimentally {g]. 
In particular, in a beautiful recent experiment j?] 'average trajectories' for photons deduced 
indirectly from the WM of momentum have been observed. In our scheme we introduce 
instead 'weak trajectories' (WT) by measuring directly the position of the wavefunction 
interacting weakly with a set of measuring apparati. We will see that in the semiclassical 
regime the only WT compatible with the positions of the pointers are the classical paths. 

Let \ip{ti)) be the initial state of a dynamical system whose evolution is governed by a 
(possibly time-dependent) Hamiltonian H{t). Let be the initial state of a measurement 
apparatus with indicating the position of the pointer. The meter wavefunction (Rfc| 0^) 
is tightly localized around some position R^,; for convenience we choose a Gaussian wave- 
function (2/7rA2)i/2e-(R-fc-R2) (we work from now on in a 2D configuration space and use 
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atomic units throughout). The local coupling between the meter and the system is assumed 
to take place during a small time interval Tk- The time- integrated interaction is taken as 
Ik = gr ■ R;j^((4A)^ ~ k ~ Rfc|^),where g is the effective coupling strength and the last term 
is a unit-step function accounting for the short range character of the interaction (this term 
will be implicit in the rest of the paper). Assume now we have a set of measurement apparati 
k = 1, ...,n each positioned around a preassigned R°. The initial state of the system and 
measuring devices |^I/(ti)) = |^(tj)) 11^=1 evolves at time t/ to 

\^{tf)) = exp -t [U {tf, t„) + J„ + ... + J2 + U{t2, h) 

n 

+h + U{hMm^))ll\<Pk) (2) 

k=l 

where U{tk+i,tk) = //j.*"^^ H{t')dt' is the self evolution of the system and tk is taken as the 
interaction time (ie Tk ^ tfc+i — tk). 

At time a standard projective measurement is made in order to postselect the system 
to a desired final state \x{^f))- Expanding each Ik in Eq. ([2]) to first order in the coupling 
g leads to 

n 

n(R.i(x(t/)i*(t/))^(x(t/)i^(t/)) 



k=l 

n 



W exp [-ig {v{tk))^ ■ Rfc] 0fc(Rfe, R^) (3) 
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where {Y{tk))y^ is the weak value [Eq. ([T])] given here by 

... ^^ _ (Xfa)|r|^(tfc)) 

with r = XX + yy. Eqs. (l3])-(jl]) indicate that as a result of the interaction that took place 



at tk, each meter wavefunction 0fc(Rfc,R°) wi 
{^itk))w -^^ standard WM scenario 



Hncur a phase-shift given by the weak value 
this phase-shift appears as a shift in the 



momentum space of each meter. Since Eq. ([3]) holds provided g and A are very small, the 
momentum space meter wavefunction will be broad, meaning a high number of events must 
be recorded in order to observe each shift. 

The structure of {r{tk))^r deserves a special comment. Each {r{tk))^r is defined a.t t = tk 
with an effective preselected state li^itk)) = U(tk,ti) I'lpiti)) being the initial state propa- 
gated forward in time and an effective postselected state (x(^fc)l = {x{'tf)\U{tf,tk) being 




FIG. 1: Time evolution of the wavefunction initially {ti = 0) given by Eq. dS]) with fq = and the 
initial mean momenta Pj,j = I, II, III taken as indicated by the arrows in panel (a). The reference 
classical trajectories /, H and M are shown resp. in black, dashed blue, and orange, (a) shows 
the wavefunction at ti = 0.7, (b) at t2 = 2, (c) at ts = 3.15 (after the wavepackets cross the origin) 
and (d) at tf = 3.65, the time at which postselection is made 8|. 



the postselected state evolved backward in time; this property illustrates the close relation 
between WM and time-symmetric formulations of quantum mechanics j4|. Note that con- 
trary to the usual definition of weak values, the effective pre and postselected states defining 
(r(tfc))^ cannot be chosen: only the initial and the final states can be freely set. A weak 
value at some intermediate time tj. reflects the interaction with the kth meter given the 
unitary evolution of the preselected and postselected states of the system. We can therefore 
envisage the set {t^, {r(tk))yy} as defining a weak trajectory of the system evolving from an 
initial state to a final postselected state as recorded by the meter devices positioned at R°, 
k = 1, ...,n. 

For an arbitrary quantum system a WT will typically reflect the space-time correlation 
between the forward evolution of the preselected state and the backward evolution of the 
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postselected state at the positions R° of the measurement apparati. Although obtaining 
this type of information is certainly of interest in general quantum systems, the notion of 
weak trajectories is particularly suited to investigate the evolution of a quantum system 
in the semiclassical regime. In this regime a typical wavefunction evolves according to the 
asymptotic form of the path integral propagator 9|, 



{2i7rfi) 



det ^ ' 



dvdv' 



1/2 



exp {iSd{r, r', t)/n- ind)} ^/'(r', 0) (5) 

where cl runs on the classical trajectories connecting r' to r in time t (from now on we 
set ti = 0) and the term between {...} is the semiclassical propagator obtained from the 
asymptotics {Sd ^ h) of the path integral form of the evolution operator U{t, 0). Sd is the 
classical action and fid the topological index of each path. Working out the full semiclassical 
propagation is often a formidable task, especially as the number of trajectories proliferate 
in the regimes where the semiclassical approximation holds. However if the initial state is 
well localized, the semiclassical propagation can be simplified by linearizing the action jl^ 
around an initial and a final reference point linked in time t by a central classical trajectory, 
the guiding trajectory. Linearization is particularly relevant if tp{r', 0) is chosen to be a 
localized Gaussian 

^ro,po(r, 0) = e-(-o)^/^^e--(--)/^ (6) 

The initial reference point is the maximum of the Gaussian, the linearized action in Eq. ([5]) 
is a quadratic form, whereas the determinant prefactor becomes a purely time-dependent 
term that can be written in terms of the stability matrix of the guiding trajectory. This 
linearized information along the central reference trajectory effectively replaces the sum over 
cl. 

The integral ([5]) can then be performed exactly: the result is of the form 
7/.ro,po(r,t) = Tr [A(t)]e-('-'iWHMW+.N(t)].(r-qW) 

gip(t)-(r-q(t))/figiS(q(i),ro,t)/fi (^j^ 

where (q(i:),p(t)) are the phase-space coordinates of the guiding classical trajectory with 
initial conditions (fq, Po) and A, M and N are time-dependent matrices depending solely on 
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the stability elements of the guiding trajectory. The advantage of working in the linearized 
regime is that by picking a preselected state of the form 

(r| HU)) = ^9V^ro,p,(r,^i) (8) 
j 

i.e. a superposition of Gaussians (|6]) launched in different directions pj, one is dealing 
conceptually with the type of problem defined by the semiclassical propagation ([5]) with a 
simplified and perfectly controlled dynamics. The evolution operator U{ti,ti) of Eq. ([2]) 
propagates each term of Eq. (jH]) along the relevant guiding trajectory yielding at ti the 
superposition of evolved states each given by Eq. J2j CjV^ro.Pj (r, ^i) (recall that q, p and 
the matrices A, M and N explicitly depend on j). 

The postselected state will also be taken to be a Gaussian of the type ([6]) localized in the 
vicinity of a chosen point at time tf 

In the linearized approximation, finding the backward propagated state is tantamount to 
obtaining the unique solution of the form ([7]) such that '?Ar(ti),p(ti)(i", t/) = Xr^,p^(r, t/): this 
gives a wavefunction centered on the time-reversed classical trajectory having boundary con- 
ditions (r/-,p/-) at t = tf. As illustrated below, suppose that the measurement apparati are 
all localized along the guiding trajectory followed by one of the branches j of the preselected 
state, so that the positions obey R° = (\j{tk). Assume for the moment that all these R° lie 
at positions where the overlap between the different branches of the system wavefunction is 
negligible. The weak values (jlj) can then be computed exactly. Two situations are then pos- 
sible, (i) If Xrf,pf{^: tk) in the vicinity of R° the meter does not move: there is no weak 
trajectory in the neighborhood of this point, (ii) Otherwise, writing the distance between the 
meter and the position at time tk of the backward evolved trajectory as = — c{f{tk), 
the weak value {r{tk))yy = (x(tfc))(^x+ {y{tk))y^y takes the form 

{x{tk))^ = [Rfe ■ x+a^-efc ■ X + 6x (Pj(4) - P/(4)) ■ x] 

+ i [9xitk)ek ■ x+/ia; (Pi(4) - P/(4)) ■ x] (10) 

and analog expressions for {y{tk))y^. a, 6, g and h are time-dependent functions whose explicit 
forms are cumbersome though straightforward to evaluate. 
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FIG. 2: (a) The postselected wavefunction, a Gaussian localized on the guiding trajectory I, is 
shown, along with the positions of the measuring devices ^^1,2,3- Only the pointer in D3 is affected, 
while P2.3 remain still: there is no WT joining the corresponding positions, (b) The 'average 
trajectories' obtained from WM of the momentum (repeated at successive positions) are plotted in 
solid black: 9 trajectories having their final positions on and near the maximum of the postselected 
state are shown, along with the reference trajectories (in faded colors). 

The structure of Eq. (fTU]) emphasizes the special role of the postselected state with 
(rj,p/) chosen such that the backward evolved trajectory simply retraces the guiding tra- 
jectory qj(t). For this special choice = and Pj(tfc) = P/(tfc) for each A;, so that each 
meter records its own position, (j{t]^)^ = R°. The WT {t^, {'r(tk))y^^} recorded by the me- 
ters thus corresponds to the guiding trajectory of the linearized Feynman propagator. For 
any other choice of postselection, {r{tk))^r can either be undefined (for those k such that 
Xr/,p/(r, tk) ~ 0) or yield a complex number depending on and the difference in momenta 
between the forward propagated preselected and backward propagated postselected states. 
The WT (where defined) becomes a list of complex numbers correlating the forward and 
backward evolution at points R°. Essentially the same situation is found when the R° are 
not placed along a guiding classical trajectory. 

For the purpose of illustration - and to avoid spurious effects due to the quality of the 
linearized approximation we will take a 2D time-dependent linear oscillator (TDLO). 
The linearized propagator for the TDLO is quantum mechanically exact, while the varying 
amplitudes capture many features of semiclassical systems with more involved dynamics. 
The Hamiltonian for the system is H = (Pj + Py)/2m + mVx(t)x^ + mVy{t)y^ where for 
definiteness we choose Vi{t) = Ki — ViCos{2uit) {i = x,y; n,v and u are constants). The 
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wavefunction ([7]) is obtained directly by employing standard path integral techniques joj; 
;he classical trajectories can be found in closed form from the solutions of Ermakov systems 
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12j . The preselected state ([H]) is taken as the superposition of 3 Gaussians at the origin 



with mean momenta as shown in Fig. 1(a). The maximum of each wavepacket then evolves 
by following the guiding trajectory, /, E or EI as shown in Fig. 1. 

Let us first set the postselected state ([9]) with r/ = q/(t/) and p/ = Piitf) and let us 
position the measuring devices Vk as shown in Fig. 2(a). The backward evolution of |x(^/)) 
simply retraces trajectory I backwards. Therefore the pointer in D3 yields according to Eq. 
fITU]) the position R3 while P2 and Pi do not move at all (no overlap with |x(^)) at any t). 
One concludes that the 'particle' went through but not through Vi and T>2- If instead 
of Vi and P2 other devices and positioned as shown in Fig. 3(a) are employed, then 
all the meters yield their respective positions as the weak outcomes: the 'particle' went 
through V^^V^ and P3. Hence one concludes (possibly by inserting additional devices) that 
the 'particle' took the WT defined by the classical trajectory /. Note that according to 
Eq. ( fTOl) there is no quantum state of the form ([9]) that can yield a WT going through "Di, 

and 1)3. This is due to the fact that there does not exist a wavepacket arriving in the 
neighborhood of at time tf that would have previously visited the neighborhoods of "Di, 

and D3. 

The last remark highlights the incompatibility between the 'weak trajectories' defined 
here by consecutive WM of the position and the 'average trajectories' (AT) defined by a 
WM of the momentum immediately postselected to a given position. By repeating these 
weak momentum measurements for different postselected positions, a velocity field is ob- 
tained. The AT are precisely the trajectories built on this velocity field. They have been 



experimentally observed recently for photons in a double slit setup 
known 



7|. It was previously 



131 ] that their dynamics is governed by the law of motion of the de Broglie-Bohm 



theory, i.e. by the probability flow, whereas the WT are generated by the semiclassical 
propagator ([5]). The mismatch [l^ between de Broglie-Bohm and classical trajectories in 
semiclassical systems hinges on the fact that when wavepackets interfere, the overall mean 
velocity field differs from the group velocity of each individual wavepacket. The mismatch 
is illustrated here in Fig. 2(b): we have computed numerically 12| the AT arriving in the 
neighborhood of = q/(tj). These AT go indeed through "Di, 1^2 and V^: starting near 
the origin, they first move in the vicinity of the guiding trajectory iZT, then travel along 
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FIG. 3: (a): The postselected wavefunction (the same as displayed in Fig. 2) is shown with the 
measuring devices T>[,'D'2 and P3. Ah these pointers indicate their own position: the evolution of 
the system along the reference trajectory / has been measured weakly, (b): Postselection now takes 
places at to = 2.84 when the wavepackets return simultaneously to the origin. The postselected 
state, chosen as a superposition of Gaussians localized around the origin with mean momenta 
pointing along the guiding trajectories is plotted along with measuring devices positioned along 
the reference trajectories /, II , M. All the meters have moved: the semiclassical sum over paths 
formulation has thus been measured weakly. 

trajectory // and thereafter 'jump' so as to move along trajectory /. 

Finally, consider choosing postselection a.t tf = to when trajectories I,E and UI first 
return to the origin, with the postselected state chosen as the superposition xo(r, to) = 
Xrj=o,po, (r, to), with Xr^=o,po, given by Eq. and poj = Pji^o), j = I, II , HI- Several 
measuring devices are positioned as shown in Fig. 3(b). By construction the backward 
evolution of xo yields a superposition of wavepackets retracing trajectories /, // and IH re- 
spectively. Therefore all the pointers will display their own position as the weak value 
indicating the 'particle' was there. This is an experimentally realizable way to catch the 
essence of the path integral approach in the semiclassical regime: weakly interacting mea- 
suring devices indicate the 'particle' takes simultaneously all the available classical paths. 
In contrast a strong projective measurement would of course yield a definite outcome on 
only one of the paths. 

To sum up, we have defined 'weak trajectories' allowing to observe the paths taken by a 
quantum system in the semiclassical regime by direct weak measurements of the position. A 
consequence worth exploring concerns the possibility of employing this scheme to reconstruct 



9 



the unknown propagator of a semiclassical system from the observed weak trajectories. The 
present setup may also be used in designing pre-post selected quantum paradoxes containing 
dynamical ingredients. 
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